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ABSTRACT 

Recent observations imply that the observed number counts of Lya emitters (LAEs) 
evolved significantly between z — 5.7 and z — 6.5. It has been suggested that this was 
due to a rapid evolution in the ionisation state, and hence transmission of the IGM 
which caused Lya flux from z = 6.5 galaxies to be more strongly suppressed. In this 
paper we consider the joint evolution of the Lya and UV luminosity functions (LFs) 
and show that the IGM transmission evolved between z = 6.5 and z = 5.7 by a factor 
1.1 < R < 1.8 (95% CL). This result is insensitive to the underlying model of the Lya 
LF (as well as cosmic variance). Using a model for IGM transmission, we find that 
the evolution of the mean IGM density through cosmic expansion alone may result in 
a value for the ratio of transmissions as high as i? = 1.3. Thus, the existing LFs do 
not provide evidence for overlap. Furthermore, the constraint i? < 1.8 suggests that 
the Universe at z = 6.5 was more than half ionised by volume, i.e. Xj^y > 0.5. 

Key words: cosmology: theory-galaxies: high redshift; statistics 



1 INTRODUCTION 

Reionization of the Intergalactic Medium (IGM) is one of 
the milestones in the history of our Universe. Following 
recombination of the primordial plasma, gas in the Uni- 
verse remained neutral until it became reionised by the first 
sources of UV-radiation. The epoch of reionization began 
with the formation of these first UV-sources, and ended 
when the entire IGM was cleared of neutral gas. The pro- 
cess of reionization was complex and is not well understood. 
The nature of the first UV-sources is unknown. Further- 
more, various feedback mechanisms, both positive and nega- 
tive, affect subsequent star, galaxy and black hole formation 
and hence the later stages of re ionization. Extensive mod- 
elling, b oth semi- analytical (e.g. ICenll2003l : IWvithe fc Loebl 
2003al lbl: iHaiman fc Holdeill2003l ) and numerical (e.g GnedinI 
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Iliev et al.l 



20061 ), 



has been performed, with the aim 
of interpreting t he few exis ting observational constraints on 
reionization (see lFan et al.ir2006a . for a review). The detec- 
tion of flux blueward of the Lya in high redshift quasars 
suggests that the Universe is fully ionised at z <.6 (e.g. 
[Fan et al. 20 j). At higher redshifts, several Sloan Digital 
Sky Sur vey (SPSS) quasa rs show complete Gunn-Peterson 
troughs (|Fan et al.ll2006a ). These troughs can be translated 
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into lower limits on the neutral fraction of hydrogen in 
the IGM, and suggest th at reionization may not have been 
completed until z m 6 jFan et ahl I2OO2I: iLidz et al.l [20o3 : 
IWvithe fc Loebll2004al : rMesinger fc Haimanll2004l2007l ). On 
the other hand, the measurement of optical depth for CMB 
photons to Thomson scattering by the Wilkins on Microwave 
Anisotropy Probe ( WMAP) of = 0.09±0.03 (|Spergel et all 
|2006| ). implies that the IGM is fully reionised out io z = 
12 ± 3 or p artially reionised ou t to even higher redshifts 
(e.g. Fig 3 of ISpergel et"alll2006l ). 

Whether reionization indeed ended at z ~ 6 is still 
an open question. By adding the derived sizes of HII 
bubbles surroundin g obse rved z — 6.5 Lya emitters, 
iMalhotra fc Rhoadd (|2006l ) have shown that at least 20 - 
50 % of the volume of t he IGM had been reionised by 2: = 
6.5. iBecker et al.l l|2006l ) suggested that the observed Gunn- 
Peterson troughs in quasar spectra at z ^6 can be explained 
without invoking an y abrupt changes in the neutral fraction 
of th e IGM (also see lMaselh et aLll2007l: iBohon fc Haehneld 
l2007h . On the other hand. iKashikawa et all (|2006l ') have 
found the Lya luminosity function (hereafter LF) to evolve 
between z = 5.7 and z = 6.5. In particular, the z = 6.5 Lya 
LF lies si gnificantly below that m easured at z — 5.7. Fur- 
thermore, [Kashi^^ra^^eF^al] (|2006l ) show that the rest frame 
UV-LF of LAEs does not evolve, within its uncertainties, 
between z = 5.7 and z = 6.5. These observations are ex- 
plained naturally, if the flux from z = 6.5 Lya emitters is 
attenuated by a larger factor than the flux from emitters at 
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z = 5.7 (lHaimanll20d3 : ISantosll20o3 ) . The increase in atten- 
uation could be interpreted as sudden change in the inter- 
galactic neutral hydrogen content, which is thought to be a 
key feature of the end of the reionization epoch. Thus, the 
Lya LF could be used to probe the last stages of reioniza- 
tion. In this paper we investigate the constraints that may 
be placed on the evolution of the IGM transmission using 
the observed Lya LPs sd, z — 5.7 and z = 6.5. As part of 
our analysis, we also investigate how the rest-frame UV-LF 
of LAEs may provide ad dition al constraints. 

Haiman & Spaan^ (|l999l ) were the first to advocate 
the use of the Lya luminosity function to constrain 
reionization, and several other pape rs have followed 
since ( e.g. iMalhotra fc Rhoa"d3 l2004l: | Le Delliou et al.l 
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CenI I2OO5I; i Le Delliou et al.1 l2006l : 
2006,; iMao et all [2OO6). The work pre- 



Haiman 
Furlanett o ct al. 
sented in this paper differs from other investigations in two 
important aspects. Firstly, our model incorporates detaile d 
calculations of the IGM transmission l|Diikstra et al.|[2007h . 
For a description of these calculations the reader is referred 
to that paper. Here, we briefly summarise the models main 
ingredients. The IGM transmission is calculated using a 
model for the gas in the IGM that accounts for clumping, 
and infall. In this model, resonant absorption of Lya 
photons by gas in the infall region (which extends out to 
several virial radii) erases a significant fraction of the Lya 
line. The model also accounts for damping wing absorption 
in cases in which a (partially) neutral IGM surrounds the 
HII bubble. Secondly, we perform statistical comparisons to 
the most up-to date observations. 

In §[2] we present our model of the Lya LF. In §[3] we 
describe the data available. The results of our comparison 
with the data and conclusions are presented in § [4] and § [5l 
respectively. The parameters for the background cosmology 
used throughout this paper are flr n = 0.24, = .76, 
Qb = 0.044, h = 0.73 and erg = 0.74 jSpergel et al.ll2006l ). 



2 THE MODEL OF THE LYa LUMINOSITY 
FUNCTION 

We use the following simple prescription to relate the Lya 
luminosity of a galaxy to the mass of its host dark mat- 
ter Mtot. The total mass of baryons within a galaxy is 
{Q,i,/Qm)Mtot, of which a fraction /* is assumed to be con- 
verted into stars over a time scale of t^ys = eoc^hub. Here, 



edc is the duty cycle and thub(z) 



3H(z) 



is the Hubble 



time at redshift z. This prescription yields a star formation 
rate of M« = f,{nt/Q,m)M/tsys- The star formation rate 
can be converted into a Lya luminosity by assuming that 
approximately two-thirds of the ionising photons absorbe d 
within the galaxy are converted into Lya l|Osterbrocklll989h . 
thus La = 0.68/i!/a,(l — /osc)Qh. Here, f^ac is the escape frac- 
tion of ionising photons, f^sc <S 1, and hva = 10.2 eV is the 
energy of a Lya photon. The total output of ionising pho- 
tons per unit mass of star formation, Qu, depends on the 
metallicity of the gas fro m which stars formed, as well a s 
their initial mass function (|Kennicuttll 19981 : ISchaere r 2003V 
In the reminder of this paper, we assume a Salpeter IMF and 



Z — 0.05^0 Q The number density of Lya emitters with Lya 
luminosities exceeding 7^ x io, is then given by 

N{>TaX La)=enc r dM-^^, (1) 

where the Lya luminosity and host halo mass. Ma are re- 
latetfl by 

^ Ma(MQ)^L 

X = 2.0 X 10^" erg s-i ' , Ta. (2) 

In this relation, Ta is the IGM transmission multiplied by 
the escape fraction of Lya photons from the galaxy. The 
function dn/dM is the P ress-Schechter (197 4) mass function 
(with the modification of lSheth et ahllioOlh . which gives the 
number density of halos of mass Al (in units of comoving 
Mpc"^). 

According to this prescription, the Lya LF has two free 
parameters {eoc, f*Ta). Thus, we exp licitly assume tha t Ta 
is independent of Ma and AT. In Dij kstra et al.l l|2007l '). we 
found Ta to decrease with increasing Ma (because infall is 
more prominent around higher mass halos) and increasing 
A'T (because a larger output of ionising photons reduces the 
impact of resonant absorption in the infall region). Con- 
versely, for a fixed eoc, increasing Ala results in an increase 
of AT, and hence an increase in Ta- Given that the range 
of observed luminosities, and hence halo masses, spans only 
~ 1 dex (§[3)1, the variation in Ta with Ma and M* can be 
ignored. In this paper we compare models for the Lya LF to 
observation as a function of these parameters and constrain 
the evolution in f*Ta between z = 5.7 and z — 6.5. 



3 THE DATA 

Detailed Lya LFs have recently been presented by Shi- 
masaku et al. (2006) for LAEs observed at z = 5.7, and 
by Taniguchi et al. (2005) and Kashikawa et al. (2006), for 
LAEs observed ai z = 6.5. Some caution must be exer- 
cised when comparing our models to this data. First, the 
observed luminosities have been derived from the observed 
fluxes by assuming that all Lya emerging from the galaxy 
was transmitted by the IGM. However, this is very un- 
likely. In fact, it has been shown that under reasonable 
model assumpti ons, the IGM transm its only a fraction of 
Ta = 0.1 - 0.3 (iDiikstra et al.1 120071 ') of Lya photons, even 



^ The uncertainty in the relation between Lya luminosity and 
M* is large, as it depends on an unknown IMF and meta llicity 
of the gas from which the galaxies formed llSchaeredliooi) . For 
the present work, because of the degeneracy between /* and the 
prefactor in Eq 2., the exact relation (and therefore the assumed 
metallicity) is not important. 

^ Although it is reasonable to assume that the star formation 
rate-and thus Lya luminosity-in a galaxy is related to the total 
amount of gas available, and thus to Mtot , no direct observational 
evidence exists to support this assumption. The fact that we ob- 
tain good fits to the data for reasonable model parameters (§0 
is promising, and suggests that our model provides a reasonable 
description. However, even if one dismisses the assumption that 
Lya luminosity relates to halo mass, then it is still possible to 
obtain constraints on the change of the IGM transmission with 
redshift using the rest-frame UV-LF of LAEs (as is discussed in 
§113. 
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Figure 1. Left Panel: Joint constraints on eoc and f*Ta- The likelihood contours at 64%, 26% and 10% of the peak likelihood are shown 
at 2 = 5.7 as dashed contours with those at z = 6.5 overlaid as solid contours. Although the likelihood peaks at different values for eoc 
and fstTa sA, z = 5.7 and z = 6.5, the contours at both redshifts fill the same region of the (edCi /*7a)-plane. A strong degeneracy exists 
between eoc ^^nd f*Ta (see text). Right Panel: Slices through the likelihood surface at eoc = 0.06 and eoc = 0.06 are shown at z = 5.7 
{black solid lines) and at z = 6.5 {red dotted lines), normalised to a peak of unity. For a given edCi f*Ta peaks at very similar values at 
both redshifts. Assuming that /» remains constant between z = 5.7 and z = 6.5, this suggests that the IGM transmission evolves only 
weakly between z = 5.7 and z = 6.5. 



when highly ionised. We therefore replace the quoted ob- 
served luminosities, Z/Lya hy TaX LLya(as in Eq[T]and Eq[2]). 
Second, the error bars on the data are uncertain. For exam- 
pl e, the open circle s show n in our Figure [2] and Figure 5 
of iKashikawa et al.l (|2006l ) denote the LF derived from the 
raw counts of their spectroscopic sample+additional photo- 
metric sample. The filled circles are corrected for detection 
incompleteness. The error-bars on the filled circles denote 
Poisson-errors. Additionally, cosmic variance is expected to 
add a varia nce of 30%(50%) at th e lowest (highest) lumi- 
nosities (e.g. ISomerville et al.]|2004l ) and enters as an overall 
shift of the Lya LF in the vertical direction. As a result, 
cosmic variance could be (partly) responsible for the off-set 
in the LFs between z = 5.7 and z = 6.5. In the first part 
of our analysis we consider the Lya LF alone, and ignore 
cosmic variance. Allowing for cosmic variance would enlarge 
the errobars on model parameters and strengthen the con- 
clusion that the evolution of the Lya LF alone, can be fully 
accounted for by the evolution of the mass function of dark 
matter halos. In the second part, we remove the effect of cos- 
mic variance by considering the joint evolution of the Lya 
and UV-LFs. 

Furthermore, the lower luminosity points (T^I/Lya ^ 
5x 10**^ erg s"^) suflter more from detection incompleteness, 
as evidenced by the larger offset between the corrected and 
uncorrected data points. In addition, the IGM is expected 
to transmit only 10-40% of the Lya flux for these galaxies, 
and the detection of these low luminosity Lya emitters re- 
quires an unusually large IG M transmission {Tg ^ 1, or un- 
usually large Lya luminosity iDiikstra et al1l2007h . It follows 
then that these galaxies must comprise only a (small) sub- 
set of the true sample. However, the exact uncertainty this 
introduces is not known. We account for this uncertainty by 
enlarging the adopted error-bars at X^Z/Lya ^ 5 x 10'*^ erg 
s^^ by factors of 4 and 3 for the lowest and second lowest 
luminosity point, respectively (on the filled circles). These 
enlargements were chosen so that the 1 — a error bars en- 



close the data points prior to the adjustment. The choice of 
error bars on these low luminosity points affects the values 
of our best fit model parameters, but not our conclusions 
regarding evolution in transmission of the IGM. 



4 RESULTS 

We calculated the Lya LF for a grid of models in 
the (edc, /*'7a)-plane, and generated likelihoods C[P] — 
exp[— 0.5x^], where = J^f^'^"'" (modeh — da.t&i)^/af, for 
each model. Here, data^ and at are the i"" data point and 
its error, and model; is the model evaluated at the i^^ lumi- 
nosity bin. The sum is over A^data ~ 6 pointsl3 

Likelihoods were determined for each redshift indepen- 
dently. In the left panel of Figure [T] we show likelihood con- 
tours in the (eoc, /*'7a)-plane at 64%, 26% and 10% of the 
peak likelihood at z = 5.7 (as dashed contours) and z = 6.5 
(overlaid as solid contours). The likelihood contours at both 
redshifts fill the same region of the (edc, /*Zi)-plane. The 
likelihoods peak at {euc , .f*Tc,) ~ (0.2,0.03) at z = 5.7 and 
at (edc, f*Ta) ~ (1.0, 0.1) at z = 6.5, with a strong degener- 
acy between euc and f,Ta. This degeneracy arises because a 
lower duty cycle corresponds to a lower tsys, which requires 
a lower /« to achieve the same star formation rate for a fixed 
value of Mtot • 

^ The data points represent a cumulative number of galaxies with 
Lya luminosities exceeding Z/Lyct- The data points are therefore 
not completely independent. However, the number of galaxies 
more than doubles between most bins, and the errorbars do not 
overlap. This indicates that although cumulative, the points are 
close to independent. 

* Reducing and ftTa by the same factor yields the same 

Lya luminosity. However, because the galaxies are visible for a 
shorter time, this still results in the prediction of fewer objects. 
This causes the deviation from a 45-degree degeneracy-line in 
the (edCi /*'^)-pl3.ne; changing cdc hy 1 order of magnitude 
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Figure 2. The Lyo luminosity function (LF) at 2 = 5.7 and 
2 = 6.5. The red squares and blue circles represent data from Shi- 
masaku et al (2006, z = 5.7) and Kashikawa et al (2006, z = 6.5). 
Overplotted as solid lines are the best fit models. For complete- 
ness, we have shown the z = 7.0 LF {black dashed line), taking 
our best fit model to the z = Q.5 data and evolving only the mass 
function between z = 7.0 an d z = 6.5. The tr iangle represents the 
2 = 7.0 galaxy discovered bv llve et al. Hiooi), which is consistent 
with the model LF. 



In Figure [2] we show the best fit models ai z — 6.5 and 
z — 5.7. Clearly, the models provide an adequate fit to tlie 
data. We also sliow tiie model z = 7.0 LF (black dot-dashed 
Ime), under the assumption that only the halo mass function 
evolved between z = 6.5 and z = 7.0. The triangle represents 
the densi t y imp lied by the z = 7.0 galaxy discovered by 
live et al] ( 20061 ) . This density is consistent with the model 
LF at the 1 — a level, suggesting that no drastic change in the 



IGM transmission is required between z — 6.5 and z = 7.0. 

At 2 = 5.7 the duty cycle was found to lie in the 
range euc = [0.03, 1.0], while at z = 6.5 this range was 
edc = [0.1, 1]. The likelihood in eoc therefore extends over 
more than 1 dex in each case. As a result, the a-posteriori 
probability for euc is quite sensitive to the choice of a prior 
probability. Thus the data do not really constrain the duty 
cycle. Therefore, in the right panel of Figure [T] we show 
slices through the likelihood surface at eoc = 0.06 and at 
eoc ~ 0.6 respectively, normalised to a peak of unity. The 
distribution of ft,Ta at z = 5.7 and z = 6.5 are shown as 
the black solid and red dotted lines. For fixed edc, the most 
likely values of f*Ta axe common at both redshifts. Assum- 
ing that /, does not evolve between z = 5.7 and z — 6.5, it 
therefore follows that the IGM transmission need not evolve 
significantly between z = 5.7 and z = 6.5 in order to explain 
the evolution of the LF. Thus the evolution in the observed 
LF may be attributed to the mass function of dark matter 
halos alone. This is explored in more detail in § 14.11 

Note that the best fit values for /.T„ = [0.01 - 0.05] 
meet our prior expectations. For /* = 0.1, we find Ta = 
0.1 — 0.5 , which corresponds t o the range of transmissions 
found in iDiikstra et al] (|2007h . In our best-fit models, the 
mass range of Lya emitters is Mtot ~ 4 — 42 x 10^° M© at 
z = 5.7 and A/tot ~ 6 - 32 x 10^° Mq at 2 = 6.5. 



is compensated for by changing f*Ta only by half an order of 
magnitude. 
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Figure 3. Likelihood contours at 64%, 26% and 10% of the peak 
likelihood in the (^,57/7i,65i /*^,57)-pl3.ne, for three different 
duty cycles Cdc- The models that provide the best fit to the 
observed LFs have 7^,57/7^,65 ~ 0.8 — 1.5. Therefore, the best 
fit-models do not favour a large change in the opacity of the IGM 
to Lya photons emitted by galaxies. 



4.1 Evolution of IGM Transmission 

In this section we investigate the evolution of the IGM trans- 
mission, Ta in more detail. We found that the observed LFs 
can be described by models over a large range in duty cycle 
EDC- Here, we fix the duty cycle and /* to be common be- 
tween 2 = 5.7 and 2 = 6.5, and vary the parameters /*'7^,57 
and R = Ta^T jTafih^ where Ta,w and Tafih are the IGM 
transmission at 2 = 5.7 and 2 = 6.5, respectively. We com- 
pute model LFs on a grid in the (f?, /t7^,57)-plane and si- 
multaneously fit to the observed LFs at 2 = 5.7 and 2 = 6.5. 

The results of this calculation are shown in Figure O 
Here, likelihood contours at 64%, 26% and 10% of the peak 
likelihood are plotted in the (i?, /«'7ij,57)-plane for 3 different 
duty cycles (eoc = 0.01,0.1 and 0.5). The likelihood-curves 
shift towards higher /*71i,57 for larger values of edc due 
to the degeneracy noted above. This figure shows that the 
ratio R = 71^,57/71^,65 is close to 1 (within <^ — a) for all 
reasonable values of eDC- Therefore, the observed evolution 
in the Lya LF can be attributed to the evolution of the mass 
function of dark matter halos. 



4.2 Constraints on Transmission from the Lya 
and UV-Luminosity Functions. 

In the previous section we derived constraints on the ratio 
7a .65 / 7(1 ,57 using a simple model for the Lya LF in which 
the star formation rate in a galaxy increases in proportion 
to Mtot, the total mass of its host dark matter halo. We 
showed that the observed evolution in the Lya LF can be 
attributed to the evolution of the mass function of dark mat- 
ter halos. As evidenced by Figure |3] our exact constraints 
are somewhat dependent on the unknown model parameter 
eDC- The goal of this section is to obtain constraints on Ta 
that are independent of any of our model parameters, and 
more generally, to obtain constraints that are independent 
of t he model that underlies the Lya LF. 

iKashikawa et al] (|2006l ) found that the rest-frame UV- 
LF of LAEs does not evolve between 2 = 5.7 and 2 = 6.5. 
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Figure 4. The likelihood distribution of the ratio Ta.65/'^a,57 
derived from the Lya + rest-frame UV luminosity function. Three 
curves are shown: the solid line shows the likelihood under the 
assumption that the UV-LFs of LAEs are identical at z = 6.5 
and 2 = 5.7. The dotted/dashed line shows the likelihood under 
the assumption that the UV-LFs of LAEs at z = 6.5 lies higher 
by a factor of 0.7 and 1.2, respectively. The grey thick region at 
R < 1.3 shows the range of ratios that can be expected if the 
Universe remained fully reionised between z = 6.5 and z = 5.7 

(§[13}. 



This appears to contradict the conclusion that the evolution 
of dark matter halos (possibly in combination with cosmic 
variance) caused the observed evolution of the Lya LF. In- 
deed, if the rest-frame UV-LF of LAEs is identical at z = 5.7 
and z — 6.5, then the entire difference between the Lya LFs 
at 2 = 5.7 and z = 6.5 must be due to a change in either 
Ta or in the ratio La/Luv (where Luv is rest-frame UV 
luminosity of a galaxy) . In other words, a lower IGM trans- 
mission (or La/Luv ratio) at z — 6.5 would shift the Lya 
luminosity function to the left relative to the z — 5.7 Lya 
LF. Assuming for the moment that La/L\jy remains con- 
stant, we quantify the ratio R = T^.es/Ta.s? that is favored 
by this constraint. While the analysis could be done for in- 
dividual galaxies, we consider the sample as a whole, by tak- 
ing the model parameters (/*,eDc) of our best-fit z — 5.7 
model. We then obtain z = 6.5 Lya LFs for a range of 
Ta,e5 by scaling the luminosities in the z — 5.7 model by 
a factor of X^.es/Ta.s?- In Figure 3] we plot the likelihood, 
C[P] = exp[— 0.5x^], as function of the ratio Tafi^lTa.w ■ If 
the UV-LF remained constant between z = 5.7 and z — 6.5, 
then we find R = lAt° l (95% CL). The preferred ratio de- 
creases if the UV-LF at z = 6.5 lies below that at z — 5.7. 
Existing data is inconclusive regarding the precise evolution 
of the UV-LF between z = 5.7 an d z = 6.5: from the size 
of the error-bars in Figure 7 of jKashikawa et all l|2006l ). we 
find that the best-fit z = 6.5 UV-LF to lie higher by a factor 
of 0.7 ± 0.2 (95%) than at z = 5.7 if all data points brighter 
than Muv = —20 are used, while it lies higher by a factor 
of 1.2 ± 0.4 (95%) if the two data points at Muv ~ -20.2 
are ignored. Figure |3] shows the the likelihood of the ratio 
'7q,65/'7q,57 assuming that the z = 6.5 LF lies higher by a 
factor of 0.7 (1.2) as the dotted (dashed) line. Assuming for 
simplicity that the ratio of the z = 6.5 and z = 5.7 UV-LFs 
has a flat likelihood distribution between 0.7 and 1.2, then 
it follows from Figure g] that 1.1 < 71 < 1.8 (~ 95% CL). 



We quantify in more detail below how this compares to the 
evolution of the transmission in a fully (§231 and partially 
(i) 14. 4p reionised universe. 

The relative normalisations of the UV & Lya LFs may 
evolve for reasons other than transmission. Either (1) Luv 
increases with redshift, or (2) La decreases with redshift. 
(1) could be caused by decreasing dust abundance towards 
higher redshift. However, the observed Lya LF requires that 
Lya is affected less by this possible evolution of the dust 
content, which in turn requires fine tuning of how the dust 
is distributed. (2) could be caused by an increasing escape 
fraction o f ionising photons with redshift as La OC (1 /"esc). 
Recently. Ilnoue et al.l l|2006l ) found that /esc is roughly con- 
sistent with the value /esc ~ 0.1 at 4 <^z <_Q. In order for 
La to in increase by a factor of ^ 1.1 between z — 6.5 and 
z = 5.7 would require /esc ^ 0.2 at z = 6.5. 



4.3 Cosmic Expansion and the Evolution of Ta 

Following the analysis presented in iDiikstra et al. I l|2007h . 
we calculate the IGM transmission (Ta) at z = 5.7 and 
z = 6.5, and the expected value of the ratio of transmis- 
sion R = Ta,57 /Ta,65- We Calculate Ta for a model in which 
Aftot = lO^^Af© (the centroid of our best-fit mass range 
for the host halos of the LAEs §|11), and A-/* = lOA/o/yr. 
The Lya line before scattering in the IGM is assumed to 
be Gaussian with a standard deviation of v^hc- The den- 
sity and veloci t y pro files of gas in the IGM are described in 
IDiikstra et al. 1 (20071 ). Note that the precise value of Ta dif- 
fers with quantities such as A/tot , Af« , and other parameters 
such as the assumed width of the Lya line, and the pecu- 
liar velocity of the ga laxies relative to the surrounding IGM 
(jPiikstra et al.l [20071 ). However, as long as these quantities 
do not vary significantly between the samples at z = 5.7 and 
z = 6.5, these uncertainties should not strongly affect the 



calculation of the ratio R. 

We consider two cases: in case I the photoionisation 
rate due to the externally generated, ionising background 
is set to 10~^^ s~^ as deri ved from quasar absorption spec- 
tra (e.g. iFan et "al]|2006bl ). In case II clustering of nearby, 
undetected, sources boosts the photoionisation rate consid- 
erably. We include the effects o f clustering using the pre- 
scription in IDiikstra et al.] (|2007l ). Clustering of sources and 
the corresponding boost in the local ionising bac kground oc- 
curs naturally in hierchical models ( [ Wvithe fe^ Locb 200f^). 
In both cases, the only difference between models at z — 5.7 
and z — 6.5 is the mean density of baryons in the uni- 
verse. Figure[5]shows the observed Lya lines after processing 
through the IGM. The horizontal axis shows the normalized 
frequencjjfl x = {v — Ua)/Auo- Note that x — corresponds 
to the true line center. The flux density on the horizontal 
axis is in arbitrary units. Here, case I (no boost) is shown in 
the left panel and case II (with boost) is shown in the right 
panel. The labels denote the values of R. 

Figure [5] shows that the infalling gas erases part of the 
Lya line redward of the Lya resonance, and produces a 



^ Here x = — Va)l I^vd^ where Aj/£) = j/ctfth/c Here Vt}^ = 
\/2kg T/mp is the thermal velocity of the hydrogen atoms in the 
gas, kg is the Boltzmann constant, T the gas temperature, nip 
the proton mass. 
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Figure 5. Lya spectra of LAEs at z = 5.7 (black line) and z = 5.7 {grey line) after processing through a fully reionised IGM. In the 
left/right panel the local ionising background is not/is enhanced due to the clustering of nearby undetected sources. The figure shows 
that R = ^,57/'?^, 65 iiia.y be as large as i? = 1.3 without invoking any drastic change in the ionisation state of the IGM. 



sharp cut-off in the Lya line at a:: > —10 fsee lPiikstra et aP 
I2OO7I . for a more detailed discussion of this). For case I, 
7? = 1.03, which is outside the acceptable range found in 
§ 14.21 For case II however, the boost in the ionising back- 
ground reduces the neutral fraction in the IGM, which allows 
a small fraction of Lya to 'leak' through the IGM. In this 
case the total fraction of Lya that is transmitted depends 
quite strongly on the mean density of hydrogen, which re- 
sults in i? = 1.3, well within the acceptable range for R and 
less than 2(t away from the best-fit model if the UV-LF of 
LAEs were identical at z = 5.7 and z = 6.5. 



4.4 Constraints from R = Ta,57 /Tafi5 < 1.8 

In the previous section we found that the ratio of transmis- 
sions R = Ta^57 /Tafis cxpcctcd for a fully reionised IGM 
out to z ^ 6.5 can be as large as R = 1.3. Following the 
analysis presented in lDiikstra et al.l (|2007f ). we calculate the 
expected value of R for cases in which galaxies at 2; = 6.5 are 
embedded in bubbles of reionised gas of radius Rmi, within 
an otherwise neutral IGM. Prior to overlap, HII bubbles 
are expected to be generated by clusters of ionising sources. 
Their sizes theref ore greatly exceed that o f HII regions of in- 
dividual galaxies ('Furl anetto et al.ll2004al lbl: IWvithe fc Loebl 
[2004b; Zahn ct al. 2007f). The associated value of 7? is shown 
in Figure |6] as a function of bubble radius for cases in which 
the local ionising background is boosted (not boosted) by 
nearby undetected sources as the dashed (solid) line. 

Figure|6]shows that R decreases with bubble size (-Rhii), 
as the damping wing absorption becomes less important. 
The minimum bubble radius shown corresponds to the case 
of an isolated galaxy in its own HII region. For large bub- 
ble sizes the ratio levels off for the boosted (not boosted) 
case at R — 1.3 (7? ~ 1.03). These are the values quoted 
in § 14.31 in reference to an ionised IGM at higher density. 
The light grey area shows the range R £ [1.1 — 1.8] that was 
preferred by the combination of the UV and Lya LF (ij |4.2|) . 
According to Figure [6l the constraint _R < 1.8 translates to 
a minimum bubble radius of -Rhii ^2 pM pc, which trans- 
lates to -Rhii <;15 cMpc. In the model of iFurlanetto et aU 
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Figure 6. The expected value of -R = TJi, 57/7^^,65 for cases 
in which galaxies at z = 6.5 are embedded within bubbles of 
reionised gas of radius -Rhii, surrounded by a fully neutral IGM. 
The grey area atl.l<-R<1.8 denotes the region of parameter- 
space that is favored by the joint UV and Lya LPs. The dashed 
(solid) line show cases in which the ionising background is (is 
not) boosted by undetected surrounding sources. Existing data 
suggests that -Rhii ^2 pMpc, which implies that the Universe is 
more than half ionised by volume at z = 6.5 (see text). 

(l2006t ) the characteristic Hll-bubble size at a given epoch is 
related to the globally averaged ionised fraction of the uni- 
verse by volume, xi^v at that epoch. If we assume that the 
LAEs at z = 6.5 reside in Hll-bubbles that are of charac- 
teristic size at z — 6.5, then the constraint -Rhii ^15 cMpc 
translates to a;i,v ^0.8. However, the actual lower limit on 
Xi^v is weaker: if the Universe were truly more than 80% 
ionised by volume at z = 6.5, then the neutral IGM sur- 
rounding the HII bubble surrounding the LAEs would be 
filled with other ionised bubbles, which would reduce the 
damping wing optical depth of the IGM, which would en- 
hance Tafis and thus reduce R. If the damping wing optical 
depth is reduced by a factor of 2 (which is roughly the case 
when the universe is half ionised by volume), then we find 
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our lower limit on bubble size to be -Rhii <;6 — 10 cMpc, 
which translates to v >0.5 — 0.6. Note that this lower 
limit on Xi^v would also be obtained directly from the con- 
straint -Rhii J; 15 cMpc, when the impact of feedback from 
clustered sources is in cluded, which modifi es the relation be- 
tween Rmi and Xi,y (jKramer et al.l[2006[ ). In summary, the 
upper limit on R suggests that a more reasonable lower limit 
on the ionised fraction of the Universe is Xi v >0.5. 



5 DISCUSSION & CONCLUSIONS 

Recent observations have shown that the observed number 
counts of Ly g emitters evolve signifi cantly between z — 5.7 
and z = 6.5 ('Kashika wa et all 120061 ). It has been suggested 
that this evolution could be due to a significant change in the 
ionisation state of the IGM during this short time interval. 
In this paper we have investigated the constraints that may 
be placed on the evolution of the IGM transmission using 
the observed UV and Lya LPs sd, z — 5.7 and z = 6.5. 

We used a simple prescription to relate the Lya lumi- 
nosity of a galaxy to the mass of its host dark matter A/tot 
and found that such a model can reproduce the data quite 
well. Using this model, we have shown that the observed Lya 
LPs at z = 5.7 and z = 6.5 are best described by a model 
in which the IGM transmission evolves only weakly between 
these two redshifts. In fact, it is possible to attribute the 
observed evolution in the Lya LP entirely to the evolution 
in the mass function of dark matter halos. The presence of 
cosmic variance in the observations strengthens this conclu- 
sion. 

However, the observed rest-frame UV-LP of Lya emit- 
ters appears not evolve between z = 5.7 and z = 6.5. Ac- 
counting for the co-evolution implies that the observed evo- 
lution of the Lya LP may indeed be due to the evolution 
of the IGM transmission between z — 6.5 and z = 5.7. 
We find in this case that the ratio of transmissions is 
1.1 < R = Tc,57/Tcfi5 < 1.8 (~ 95% confidence levels). 
This result is insensitive to the underlying model of the Lya 
LP (as well as cosmic variance). However, we find that the 
ratio of transmissions expected for a fully reionised IGM 
out to « > 6.5 can be as large aa R — 1.3. Thus as with 
consideration of the Lya LPs alone, the observed evolution 
is consistent with no change in the ionisation state of the 
IGM. The existing LPs therefore do not provide evidence 
for overlap between z = 5.7 and z = 6.5. Purthermore, the 
upper limit R = Xit.sT/Ta.es < 1.8 implies that the Uni- 
verse at z = 6.5 was more than half ionised by volume, i.e. 
Xi^v > 0.5. 

The present paper is based on observed LPs that were 
derived from spectroscopic observations of 17 z = 6.5 and 
28 z = 5.7 Lya emitters. A larger sample of Lya emitters 
at both redshifts will allow more stringent constraints to be 
placed on simple models such as those presented here. In 
particular, the existing LP covers 1 dex in luminosity only. 
A larger range in luminosities will be useful for (i) breaking 
model degeneracies and determining the duty cycle and (ii) 
determining the evolution of the IGM transmission to 
greater accuracy, especially when Lya LPs are used in com- 
bination with the UV-LP of Lya emitters. This will provide 
constraints on the epoch of reionization which are indepen- 
dent of those derived from quasar absorption studies. An 



improved determination of the observed Lya LP will also 
allow constraints to be placed on more sophisticated mod- 
els, that account for the impact of the i nterstellar medium 
on th e Lya emission from galaxies (as in lHaiman fc Spaanj 
Il999l ). include galaxy clustering and account for scatter in 
quantities such as f*Ta, euc etc. These issues will be ad- 
dressed in future work. 
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